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Abstract 

We construct a family of quasi-solvable quantum many-body systems by an algebraic method. 
The models contain up to two-body interactions and have permutation symmetry. We classify 
these models under the consideration of invariance property. It turns out that this family includes 
the rational, hyperbolic (trigonometric) and elliptic Inozemtsev models as particular cases. 
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I. INTRODUCTION 



New findings of solvable or integrable models have stimulated development of new and 
wide research directions and ideas in both physics and mathematics. The discovery of 
quasi- solvability in quantum mechanics is a typical example. By quasi-solvability we 
mean that a part of the spectra can be solved, at least, algebraically 1 . One of the most 
successful approach to construct a quasi-solvable model is the algebraic method introduced 
by Turbiner in 1988 |2| , in which a family of quasi-solvable one-body models was constructed 
by the sl(2) generators on a polynomial space. This family was later completely classified 
under the consideration of the GL(2, M.) invariance of the models 0,0]- Recently, this famih 
r e been paid much ; 

ill, El EI El EI El- Several attempts were made to construct quasi-solvable many-body 



have been paid much attention to in the context of A/"-fold supersymmetry 



sntly this family 



models by naive extension to higher-rank algebras. Especially, construction of two-b ody 
problems by the rank 2 algebras was extensively investigated jl EI EI EI HI 21, H- 



These approaches however led to Schrodinger operators in curved space in general and could 
hardly apply to M-body (M > 2) problems. 

In 1995, a significant progress was made in Ref. (23|, where the exact solvability of the 
rational and trigonometric A type Calogero-Sutherland (CS) models (H HI HfjJ for any 
finite number of particles were shown by a similar algebraic method. The key ingredient 
is the introduction of the elementary symmetric polynomials which reflect the permutation 
symmetry of the original models. The algebra for the M-body system is sl(M + 1). This idea 
was further applied to show the exact solvability of the rational and trigonometric A and BC 
type CS models and their supersymmetric generalizations and to show the quasi-exact 
solvability of various deformed CS models [H Hj|. Therefore, one can say the approach 
starting from Ref. |23| is, up to now, the most successful in investigating quasi-solvable 
quantum many-body problems. However, one has not yet known all the models that can be 
obtained by this approach. In other words, we have not obtained the classification of these 
s[(M+l) M-body models like that of the sl(2) one-body models. Recently, this classification 
problem was partly accessed in Ref. 3(| though, as was stressed by the authors themselves, 



the results depend on the specific ansatz and thus are incomplete. In this Letter, we will 
show the complete classification of the quantum many-body systems with up to two-body 
interactions which can be constructed by the sl(M + 1) method. 



II. CONSTRUCTION OF THE MODELS 

Consider an M-body quantum Hamiltonian, 

1 M d 2 

H * = ~2^^ + V(q u ...,q M ), (1) 

i=l 1 

which possesses permutation symmetry, that is, 

V( . . . , ft, . . . , qj, . . .) = V( . . . , q j: . . . , q { , . . .), (2) 



1 The term quasi-exact solvability has been widely used in this meaning. However, we keep it to express 
the case where the state vectors corresponding to the solvable spectra are normalizable. Importance of 
this distinction is explained in Refs. |12lll5| 
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for Vi 7^ j. To algebraize the Hamiltonian ((TJ, we will proceed the following three steps. At 
first, we make a gauge transformation on the Hamiltonian (Jlj): 

H N = e^H N e~^\ (3) 

The function W(g) is to be determined later and plays the role of the superpotential when 
the system Eq. (JIJ is supersymmetric. As in Eq. (|3*|). we will hereafter attach tildes to 
both operators and vector spaces to indicate that they are quantities gauge-transformed 
from the original ones. In the next, we change the variables to hi by a function h of 
a single variable; hi = h{qj). Note that the way of changing of the variables preserves 
the permutation symmetry. The third step is the introduction of elementary symmetric 
polynomials of hi defined by, 

a k (h)= K---K (k=l,...,M), (4) 

ii<— <ik 

from which we further change the variables to <7j. Then, we choose a set of components of 
the A/"-fold supercharges in terms of the above variables Oi as follows: 

p N=*Z (l<ii<"-<W<M), (5) 



where {i} is an abbreviation of the set . . . Using these A/"-fold supercharges, we 



define the vector space Vat = Hm ker pffi , which now becomes, 



V N = span {a™ 1 • ■ • <f : < E,=i^ < A/" - 1 j . (6) 

For given M and A/", the dimension of the vector space © becomes, 

^ (n + M-l)! (A- + M-1)! 
dim ^=E w!(M -l)! = (A--1)!M! • (7) 



n=0 



We will construct the system Q to be quasi- solvable so that the solvable subspace is given by 
Jus t Eq. (jfij). This can be achieved by imposing the following quasi-solvability condition [l^, 

MM, 

P$H N V N = Q for V{z}- (8) 

The general solution of Eq. (JHJ) can be obtained in completely the same way as shown 
in Refs. (Hill. As in the case of the one-body models, it is sufficient to find differential 



operators up to the second-order as solutions for since we are constructing a Schrodinger 
operator in the original variables It turns out that the general solution which contains 
up to the second derivatives takes the following form, 

M M 

Hj^f = — 2J A K \^ V E ' K \E ' M „ + 2J B K \E K \ — C, (9) 

K,\,/1,U=0 K,X=0 
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where A K \ j/JlV , B K \, C are arbitrary constants, and E K \ are the first-order differential opera- 
tors which constitute the Lie algebra sl(M + 1): 

d d 
E 0i = q^-, E i:j = (10a) 

M a 

a 



E t0 = a t E 00 = ^1 M - 1 - ^ a k — 1 . (10b) 



k=\ 



If we explicitly express the general solution (jHJ) in terms of crj, we obtain the following 
expression, 



M 



Hjv = - [A (a)a k ai - A k (a)a t + A kt (a 



d 



.2 



k,l=l 
M 



da k dai 



+ J2[Bo(v)<rk ~ B k {a)] A - C(a), (11) 



where A K , A k i, B K and C are second-degree polynomials of several variables. 

One of the most difficult problems one would come across in the algebraic approach to 
the quasi-solvable quantum many-body systems is to solve the canonical-form condition: 



i M pa 

H N = e-WV« = -\Y,h + (12) 



2 ^ dq: 

i=i 

If the Hamiltonian (fTTj) is gauge-transformed back to the original one, it in general does 
not take the canonical form of the Schrodinger operator like Eq. (JTJ) and one can hardly 
solve, for arbitrary M, the conditions under which a gauge-transform of Eq. (fTTj) could 
be cast in the Schrodinger form. This difficulty can, however, be partly overcome by the 
following observation. Suppose we can solve the canonical-form condition for an M and 
obtain a quasi-solvable M-body Hamiltonian constructed from the sl{M + 1) generators, 
which would have the following form: 

M q2 M M 





i<j 



H " = ~\ Yl Qtf + S + S V ^ *) + •■■ + v m(.Qi, ■■■,Qm). (13) 



Then, we can get a quasi-solvable M-body model with up to the two-body interactions if 
we turn off all the coupling constants of the interactions except for the one- and two-body 
ones: 

M q 2 M M 



d 2 

? + 

i=l i<j 



The resultant model (|14|) should be, when we put M — 2, identical with one of the two-body 
models constructed from the sl(3) generators. This comes from the fact that the gauged 
Hamiltonian (fTTj) constructed from the sl(M + 1) generators reduces to the one constructed 
from the st(3) generators if we put M = 2 and hi = for i > 2. 
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Therefore, as far as up to two-body interactions are concerned, it is sufficient to solve the 
M = 2 case by virtue of the permutation symmetric construction. We have found that we 
can actually solve the canonical-form condition for M = 2 and that Hj^ for M > 2 must 
have the following expression in terms of the variables h i} 



M 



8=1 

M 



-E 

i=l 



Q{K) - 



7V-2 + (M- l)c 



_d_ 



M 



P(hj) d 
hi — hi dhi 



C(a(h)), 



(15) 



where C is given by, 

C(a(h)) 



12 
J\f 



[M - 2 + 2(M 
-1 - 



Ai 



-l^^P"^) 

i=i 

9 



i=l 



ft,,- — /i, 



+ R. 



(16) 



The P and Q in Eqs. fjTHjl and (fTT)|) are a fourth- and a second-degree polynomial, respec- 
tively: 



P(h) 
Q(h) 



a^h + CZ3/Z 3 + a<ih 2 + ai/t + ao, 
b 2 h 2 + b 1 h + b . 



(17a) 
(17b) 



Thus, there are 10 parameters a„, & n , c, P, which characterize the quasi-solvable Hamiltonian 
()15|). One can prove the quasi-solvability of the operator ()15|) by the convertibility of it into 
the form (jllj) . The function /i(g), which determines the change of variables, is given by a 
solution of the differential equation, 



h'{ q f = 2P{h{q)). 



(18) 



One may notice that the resultant Eqs. ()15j) - (fT8~j) have resemblance to those of the one-body 
quasi-solvable models constructed from sl(2) generators @,IEIH|> or equivalently, the type A 



Af-iold supersymmetric models 0, 3 JiJ ■ Indeed, we can easily see that the above results 
reduce to the one-body st(2) quasi- solvable and type A A/"-fold supersymmetric models if we 
set M = 1, where the double summation is understood as = 0. This is consistent with 
the fact that in the case of M = 1 the above procedure is essentially equivalent to that in 
the s[(2) construction of type A A/"-fold supersymmetry (I3.E^|. Under the above conditions 
(|15 )) -([18 p satisfied, the original Hamiltonian becomes the following Schrodinger type, 



H 



i=i 



M 



i=l 



( dW{q) \ 2 d 2 W{q) 



V 9qi J 



dq\ 



C{a(h)) : 



(19) 
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and the superpotential W(g) is given by, 



w*) = -£ + * ~ 1 T - 1)c f> w - °f> i* - M • (20) 

i=l ^ ^ l ' i=l «<j 

It is evident by the construction that the solvable wave functions ip(q) of the Hamiltonian 
(fTTJj) take the following form, 

^)=#?)e- W(9) , (21) 

The Hamiltonian (|19|) with Eqs. and (J2(jp is the most general quasi-solvable many-body 
systems with two-body interactions which can be constructed from the sl(M + 1) generators 

(HOD- 

Before investigating what kind of particular models emerges from the general Hamiltonian 
(|19|). we will refer to an interesting feature of the result. If the algebraic Hamiltonian 
does not contain any raising operator Eio, it preserves the vector s pac e I V for arbitrary M 
and becomes not only a quasi-solvable but also a solvable model |2lL l22j ]. In this case, it 
turns out that C(er) = C, one of the constants involved in Eq. Q, and thus the original 
Hamiltonian (J19)) becomes supersymmetric 0,H2|. A system is always quasi-solvable if it 
is supersymmetric, since the ground state is always solvable. From the above result, we can 
conclude that a system is always supersymmetric if it is solvable and all its states have the 
form (j2U). 



III. CLASSIFICATION OF THE MODELS 

It was shown that the one-body sl(2) quasi-solvable models can be classified using the 
shape invariance of the Hamiltonian under the action of GL(2,M.) of linear fractional trans- 
formations 0, 0|. We can see that the many-body Hamiltonian (jl5j) also has the same 
property of shape invariance. The linear fractional transformation of hi is introduced by, 

hi^k = ^^j (A = a*-/3 7 ^0). (22) 

Then, it turns out that the Hamiltonian (JT5j) is shape invariant under the following trans- 
formation induced by Eq. (|22|). 

M M 

H N {h) » H N {h) = + H N {h) JJ( 7 ^ + 8)-^- l \ (23) 

i=l i=l 

where the polynomials P(h) and Q(h) in the Hj^{h) are transformed according to, 

P{h) i — > P(h) = A~ 2 {-fh + 6)*P(h), (24a) 
Q(h)^Q(h)=A- l ( 1 h + 5) 2 Q(h). (24b) 

For a given P(h), the function h(q) is determined by Eq. (jl8j) and a particular model is 
obtained by substituting this h(q) for Eqs. (fTEJ). (fT^j) and (j2*Uj) . Under the transformation 
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(j24aj) of GL(2,M), every real quartic polynomial P(h) is equivalent to one of the following 
eight forms: 

1). ~ 2). 2/i, 3). 2uh\ A).2u{h 2 -l), 

5). 2^ 2 + l), 6). ^(h 2 + l) 2 , 

7). 8). U -(h* + l) [(l-k^ + l] . 

where i/, k, g-i and (73 are all real numbers satisfying 0<A;<1 and g\ — 27g 2 ^ 0. 

Thus, the quasi-solvable models (fTSjl can be classified into the above eight cases. 

Case 1). h(q) = q: 

This leads to the rational A type Inozemtsev model j^, Q, |35| . Inozemtsev models are 
known as a family of deformed CS models which preserve the classical integrability. The 
main difference between quantum and classical case is that the quantum quasi-solvability 
holds only for quantized values of the parameter, say, for integer Af, while the classical 
integrability holds for continuous values. This is one of the common features that the 
quantum quasi-solvable models share. 

Case 2). h(q) = q 2 : 

This leads to the rational BC type Inozemtsev model. The quasi-exactly solvable model 
reported in Ref. j2^| is just this case. 

Case 3). h(q) = e 2 ^: 

This leads to the hyperbolic [y > 0) and trigonometric [y < 0) A type Inozemtsev model. 

Case 4). h(q) = cosh2y / t/g: 

This leads to the hyperbolic {y > 0) and trigonometric [y < 0) BC type Inozemtsev 
model. The quasi-solvability of the special cases of the above four were recently shown in 
Ref. [36] by an ansatz method. 

Case 5). h(q) = smh2y/uq: 

This leads to a hyperbolic model being neither the Inozemtsev nor the Olshanetsky- 
Perelomov type [37|. The paper Ref. covers most of the above five models. 

Case 6). h(q) = tenay/vq: 

This leads to a trigonometric model being neither the Inozemtsev nor the Olshanetsky- 
Perelomov type. 



Case 7). h(q) = p(q; g 2 , g 3 )- 

This case includes the elliptic BC type Inozemtsev model and the twisted CS mod- 
els I39L l4fl| . The elliptic model in Ref. |4l| may be also included in this case. 

Case 8). h(q) = sn (y/uq\k)/ cn (\fvq\k): 

This leads to an elliptic model being neither the Inozemtsev nor the Olshanetsky- 
Perelomov type. 
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The two-body potentials in all the cases have singularities at g« = qj (i ^ j). Thus, each 
of the models is naturally defined on a Weyl chamber if the potential is non-periodic or on 
a Weyl alcove if the potential is periodic [37j • Cases 1-5 with v > and Case 6 with v < 
correspond to the former while the others to the latter. In the latter case, a system can 
be quasi-exactly solvable unless a pole of a one-body potential in the system exists and is 
in the Weyl alcove. On the other hand, quasi-exact solvability in the former case depends 
mainly on the asymptotic behavior of Eq. (|21j) at \qi\ — > oo. Since this behavior is in general 
not dominated by the two-body term in the r.h.s. of Eq. (|20p. most of the results on the 
normalizability of the one-body sl(2) quasi-solvable models in Ref. may also hold for our 
models. 

More details on the results presented here and further development will be reported in 
the near future [42] . 
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